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Abstract 

Cold dark matter (DM) scenario may be cured of several problems by involving self-interaction of 
dark matter. Viability of the models of long-range interacting DM crucially depends on the effectiveness 
of recombination of the DM particles, making thereby their interaction short-range. Usually in numeric 
calculations, recombination is described by cross section obtained on a feasible quantum level. However 
in a wide range of parameter values, a classical treatment, where the particles are bound due to dipole 
radiation, is applicable. The cross sections, obtained in both approaches, are very different and lead to 
diverse consequences. Classical cross section has a steeper dependence on relative velocity, what leads to 
the fact that, after decoupling of DM particles from thermal background of “dark photons” (carriers of 
DM long-range interaction), recombination process does not “freeze out”, diminishing gradually density 
of unbound DM particles. Our simplified estimates show, that at the taken parameter values (the mass of 
DM particle is 100 GeV, interaction constant is 100 -1 , and quite natural assumptions on initial conditions, 
from which the result is very weakly dependent) the difference in residual density reaches about 5 orders 
of magnitude on pre-galactic stage. This estimate takes into account thermal effects induced by dipole 
radiation and recombination, which resulted in the increase of both temperature and density of DM 
particles by a half order of magnitude. 


The models of self-interacting dark matter (DM) have aroused a lot of interest in the last time [Til. 
DM with long-range interaction (referring hereafter as to y-interaction) seems to be able to escape several 
problems of ordinary cold dark matter (CDM) scenario, such as an overproduction of subhalos and cuspy 
density profile in them [6j8, 12[|l3]. At the same time, an ellipticity of big halos is not spoiled at some 
model parameters j'8]. An enhancement of annihilation signal in the Galaxy (so called Sommerfeld-Gamov- 


Sakharov enhancement |14f[l6] ), considered for the first time (to our knowledge) in 17 , is one more possible 
bonus of the models of question. Analysis of recent observations of forming galactic cluster Abell 3827 also 

Origin of supermassive black holes can be connected with an existence 
Generally, models with dissipative form of DM as 


favours self-interacting DM 18 


of DM component with strong self-interaction 19 


sub-component find more applications 20,21 


Essential feature of cosmological evolution of y-interacting DM is a formation of atomic-like bound states 
by DM particles with opposite y-charges. If oppositely y-charged particles are particle (a) and anti-particle 
(b = a), then they annihilate, what may drastically affect their residual density 15,12,221. If the bound 
particles are different species (a and b / a) so bound state is stable, then depending on relative amount of 
bound and unbound particles, as it is obtained by the period of large scale structure formation, DM dynamics 
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is very different and whether or not the model gets constraint from observations J3 12,131. Description of 
recombination process plays a clue role here. Usually a quantum approach is used for it. However a classical 
approach, which was used, in particular, for magnetic monopoles and heavy neutrinos 122 23 


seems to 

be valid in a broad interval of parameter values. It leads to a result very different from that obtained on 
quantum level, which in the commonly accepted form does not come to the classical limit. 

Classical recombination cross section is obtained from condition that the scattered particles lose, due to 
dipole radiation, sufficient energy to get bound |24,25 , and is given by 


CTree = TrPmaxW = (4vr) 2/5 7r 


1 


p 2 V 


14/5 ’ 


( 1 ) 


where p max is the maximal impact parameter at which a pair is bound, a y is the constant of y-interaction, 
A 4 = is the reduced mass of the pair of the scattered particles with m a and my, being their masses 

(m a < m b ), v is their initial relative velocity. This cross section has a steeper dependence on velocity with 
respect to that of usually accepted quantum recombination cross section. The latter in form of Kramers’ 
formula [26| summed over all quantum levels is (valid for v <C a y ): 
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(2) 


Steeper behaviour of classical cross section leads to the fact that recombination process does not freeze 
out on both radiation dominated (RD) and (even faster expanding) matter dominated (MD) stages, and 
relative number of unbound y-charged particles falls down gradually with time. 

Classical formula ([I]) is assumed to be valid when |25] 

v < a b J 2 . (3) 


Under this condition, binding is found to occur predominantly due to multiple soft photon emission, what 
allows classical treatment. However considering on quantum level, only one-, two-photon final states are 
usually taken into account. Eq.Q can be formally deduced from condition that binding of two particles 
(i.e. when initial kinetic energy of relative motion, E re \ = pu 2 /2, is lost) happens at distance (Rb) much 
greater than the radius of the respective ground bound state (ob). In this case, action of the system, as will 
be shown, becomes much greater 1 (in units h = 1) and thus the classical approach is reasonable. 

Below we shortly discuss classical approach implications in recombining DM cosmological evolution. It 
includes estimations, some of which may seem to be simple, to trace explicit dependence of result on the 
parameters. 


Our results will basically relate to parameter region lying around fiducial values m a ~ p = 100 GeV and 
a y = 1/100 used for numerical estimations. Also for definiteness we assume (as in case of heavy neutrino 
model |22||) that before a direct annihilation of a and b (where either b = a or b ^ a), happening when 
the temperature becomes below their mass (of the lightest), T = T* ~ m a /10, y-plasma (consisting of a, 
b with their antiparticles and y-photons - massless y-force carriers) has the same temperature as ordinary 
matter. Right after annihilation, y-photon background decouples from a and b as well as from ordinary 
matter (O), while the opposite is not true. That is the y-background is no longer influenced by a and b, 
and by O through O-a, b possible coupling, but a and b are influenced by y. Starting from this moment, 
the temperature of y-background ( T y ), as of a closed system, changes as inverse scale factor, whereas that 
of ordinary photons feels also entropy re-distribution between ordinary matter components. So for y- and 
O- photons’ temperature relation one has 


Ty = K 1/3 T, 


k(T) 
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where g s ,o(T ) is the effective number of ordinary matter species (excluding y) contributing into entropy 
density (s). For the chosen numeric values, contribution of y in density at nucleosynthesis (BBN) makes 
up At 4 / 3 (T 1 MeV) ~ 0.06 from that of O-photons, what has no effect on BBN data. The temperature of 

a and 6 is equal to T y until they decouple from y, influence of O-matter (being determined by some weak 
scale interaction) is negligible (see Appendix A). After decoupling it evolves as (see Appendix A) 


Ta ~ T 2 /T c 
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Here g e , being effective number of matter species (including y ) contributing into energy density, as well as 
At are taken at T = 0.2 MeV. 

Evolution of number density of unbound a- and b- species can be approximately described by equation 
system (see Appendix B) 


irp ( (frecV ) Tjrp 

de_ f ay / ( 1 2 \ \rs (6) 

dT ~ T 2\ 3 E P™ 9 E "*j VrecV J ^ 


Here r is the number density conventionally expressed in units of s, the brackets “( )” mean averaging over 
velocity distribution of a and 6, 9 shows deviation of T a from Eq.([5]), E pa j r = E a + is the energy lost by 
a-,6-gas (thermal bath) in the result of the pair binding (that is the case of annihilation, i.e. when 6 = a, or 
when the bound systems are thermally decoupled from a-, 6-gas). The second equation takes into account 
thermal effects, caused by scattering of particles: (*) presumably slower particles to be bound “go out” 
(annihilate or decouple) of a-,6-plasma, effectively heating it; (**) scattered but unbound pairs experience 
dipole energy losses cooling plasma. Evolution is considered in terms of O-photon temperature (T). 

Here we do not take into account some recombination process details which are more appropriate for 
quantum case (such as recombination into different level bound states, red-shifting of recombination pho¬ 
tons), as well as inverse processes (which are not important for the big parameter space of question), and 
also in quantum case the second equation of ([6]) is omitted (9 = 1). 

The bound states start to form when T a becomes much lower than ionization potential / = ya y /2, 
T a < T arec ~ 1/ 10. For the chosen values, it is T arec = 0.5 MeV, at which T a = T y , what corresponds to 
O-matter temperature T = T rec = K~ l ^T aiec ss 1 MeV. Depending on behaviour of (a rec v) with T and 
epoch, recombination process should flow in different regimes. Basically, it either damps (freezes out) and 
it is initial moment (when it starts) what predetermines the residual density of free a and 6, or it “burns” 
continually with a self-adjusted rate and final moment, at which we need to know the density, defines its 
value. For classical cross section at T ~ T rec recombination process gives effect and temporarily freezes 
out, but after a-y decoupling it is restored and goes with a steady rate until the galactic stag^] Herewith, 
recombination rate r rec = n ( a rec v ) turns out to be of the same order as Hubble rate ( H ), if thermal effects 
(second equation of system ([6])) are ignored (9 = 1), otherwise (6 / 1) the ratio T iec /H slowly goes down 
with time. So the value of residual density is not sensitive to the initial moment and all early history, but 
fully defined by the final one, and is just weakly dependent on initial abundance for 9^1 (see below). 

Classical approach is assumed to be valid when Rf, as well as energy loss length-scal^J Li oss , are much 
greater than dark atomic size as- From other side, they must be much less than spacing between y-charged 
particles, L sp (a fortiori a-a(b),y interaction lengths). Rb and Ti oss can be found by taking mechanical 

Mhe moment T a ~ 7/10 can come after a-y decoupling at some parameter values and the first freezing stage is absent. 

J It is defined as the length on which 90% of initial energy is lost due to dipole radiation before the particles are bound. 
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Figure 1: The ratio Rb/as as function of impact parameter for two redshift values (left) and of temperature 
for two values of y-interaction constant (right). 


energy (kinetic plus potential in cms) loss rate equal to dipole radiation intensity 


dE 1 "n 

— - — = dA 

dt 67T 

Dipole moment can be expressed with the help of Newton’s law, d = e y r, |r| = yf j- t (braking due to radiation 


is negligible here), and then through relation dt = dr/ 
from distance between the particles r 
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Here M is the angular momentum, which can be assumed to be conserved. In the region of interest, solution 
of Eq. ([T]) can be simplified by neglecting U E” in the square root. 

Condition Rb, L\ oss <C L sp is found to be true by a wide margin for the most of parameter values at any 
epoch of question 27 . The values Rb and Li oss are of the same order of magnitude and both depend on 


impact parameter p and on energy, which falls down with the Universe expansion. Dependences of Rb/as 
from p/pmax and from the temperature of O-photon are shown in Figs (1(a), 1(b)). One can see that Rb/a-B 


at the most of p values is close to that at p = 0, therefore Rb/as{p = 0) is used in all further calculations. 


1 - ( ^ 
i r 


5/2' 


In this case Rb ~ and E ~ EE. 

Since the late period is determinative for residual density value in case of classical cross section, we 
single out late expansion moment, when the galaxies start to form, z ~ 10 (T ~ 30 K). As one can see 
from Fig (1(b)), Rb/as S> 1 there. This ratio depends on parameters as Rb/ciB oc a y p 2 ^ 5 /T^ 5 ~(after a-y 

a^/’/i/T 4 / 5 (for 0 = 1). At the chosen parameter values, classical approximation does not 
200 eV, but, importantly, it does in the late period. 


decoupling) r 
work for T > T q _ c 

Solution of the system ([6]) is given by Eqs.(25)-(27) of Appendix B. In the Fig.(2(a)) the density of 
unbound y-charged particles, obtained in different approximations, is shown as function of T. In classical 
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Figure 2: Left: relative density of unbound y-charged DM particles as a function of T, obtained with either 
classical (dashed curve) or quantum (dashed-dot curve) recombination cross sections (<7 rec ). The curve “Q- 
C” is obtained applying before and after the moment T ~ 200 eV quantum and classical a rec respectively. 
On the small plot inside, a quantum case is shown with the scaled vertical axes. Right: the region of values 
/r and a y , where classical cross section should be used for estimation of residual density of (un)bound DM 
particles. 


case, the mentioned above features are seen: sharp recombination effect freezes out around T ~ 1 MeV, 
continual recombination takes place from T ~ 0.1 MeV until the late epoch, changing slope at RD—)-MD 
transition. (All sharp breaks in the curve would be smoothed if estimated more accurately.) Quantum 
formula gives very small effect at the chosen parameters: there are unessential effect at T ~ 1 MeV and 
a slow logarithmic decline of r(T) from y-a decoupling to RD—»MD transition. The curve “Q-C” in the 
figure is obtained with application of Eq.Q in interval T > T q _ c and Eq.Q at T < T q _ c . As one can see, 
in the late period this solution comes to that obtained in pure classical approximation. A small deviation 
is caused by thermal effects (described by second equation of the system ([b])) which make the result being 
a little sensitive to the initial density (see Eq.(27) of Appendix B). In calculations the following values had 
been used: a y = 1/100, m a = 100 GeV, mb = 1 TeV, r(T rec ) = pcdm/ttW s(T = 2.7K) ~ 4.6 x 10 -13 
with pcdm ~ 1.4 keV/cm 3 being the modern CDM density. Also, quantum recombination rate was taken 
from [l], where <r rec is a little corrected as compared to Eq. ©• 

Note that the mentioned above thermal effects in classical case are found to be weak ( 6(T ) rsj T -1/23 ), 
but nonetheless noticeable. At T ~ 30 K for y = 100 GeV and a y = 1/100, 6 ~ 4. The ratio of r(T ~ 30), 
obtained with and without account for thermal effects (using only classical cross section), is ~ 3. 

Condition © for period z ~ 10 allows to outline the region of parameters values {y and a y ) when 
the classical approximation for a rec is applicable. It is shown in the Fig.(2(b)). There can exist a region 
of parameter space where calculations should be re-considered, since interaction radius ~ a y /T a exceeds 
spacing between DM particles [27!. 


Finally, we show arguments in favour of condition Eq.Q on the base of action in the Appendix C. 


To conclude, we have shown that classical recombination cross section Eq.Q can be applicable in a 
broad parameter region of self-interacting DM models, and it leads to strongly different result comparing to 
obtained with usually accepted quantum cross section. It may change (extend) the parameter region of the 
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corresponding models’ viability. 
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Appendix A 

Here we trace shortly the temperature evolution of y-charged particles, a ( b ), before recombination starts. 
They are assumed to experience energy exchange with y-background and, possibly, ordinary matter. In 
calculations we suppose that a interacts with O as Dirac heavy neutrino. 

To find T a one can formally use the first law of thermodynamics, which for a-particles can be reduced to 

3 dT 

2= ( ^ aV ) ay % + ( ^ aV )ao n °~ 3T * H ’ ( 8 ) 

where o = e,v,p,n,... are available species of O-matter, ( AEav ) . is the kinematically averaged energy 
transfer in ^-interaction multiplied by respective cross section and relative velocity, averaged in thermal 
distribution, n* is the respective number density (i = y or o). 

In all scattering processes of interest a is non-relativistic, i is ultra-relativistic (p, n turn out to be too 
suppressed in density, so play no role in heat transfer with a). We adopt Boltzmann approximation for all 
species distributions. The ay-scattering is well described by Thomson cross section over the great part of 
parameter space. Then calculation of ( AEav with accuracy ~ T/m a gives 


( AEav ) 

\ /ay 
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where ui and E are the energies of y- and a-particles respectively. 

To estimate ^ AEav ^ , assuming that a is heavy neutrino, one needs to take the cross sections of 
relevant processes ( au , aP, ae ~, ae + ), which in our limit are undistinguishable for particle and antiparticle 
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Here wi a b is the energy of incident z-particle in the reference frame where a is in the rest, Gf is the Fermi 
constant, £ e = 1 — 4£ + 8£ 2 ~ 0.50 with ^ = sin 2 9w being the weak mixing parameter. Unlike ay-scattering, 
ae-, av- cross sections depend on energy, what accounts for higher power of the temperature in the final 
expression 
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T u = T e = T was put. Coefficients T ay ~ 0.1 MeV, T ao ~ 10 MeV. Second term in the r.h. of Eq.(12) has 
no effect on solution. Excluding it, solution can be expressed in the form 


T a (T) = 


^k^T 2 


exp(T 2 /T 2 )(l-erf(T/T iVy )). 
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At T <C T ay (i.e. after a-y decoupling) 
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Appendix B 

Evolution of abundance of free a-, b- particles can be described by Boltzmann equation 

^ = -n 2 (CT rec u) - 3 Hn, 


(15) 


which is easily reduced to the first equation of the system Eq. 
r^T 3 , _ dt= i_dT (for 

one can get 


9e,: 


with the help of replacements: n = rs = 
const). For calculation one parametrizes <7 rec = cfq/v^. For (cr iec v) then 
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Since recombination rate is strongly temperature dependent (especially in classical case), thermal ef¬ 
fects of recombination process itself can be important, correcting temperature evolution. These effects are 
relevant, obviously, after a-y decoupling. 

One can take again the first law of thermodynamics, dQ = 5A + dU. One has the total particle number 
in some volume N a b = n a bV with n a b = n a + rib = 2n being their number density, the pressure p = n a bT a - 
Expansion of the Universe is treated as a work of gas: 5A = pdV = n a bT a 2>HVdt. Inner energy gain is 
dU = %d(N a bT a ) = ^N a bdT a + ^T a dN a b. Here we assume that dN a b = —2 (cr rec u) n 2 Vdt as if the recombined 
pairs disappear from o-, b- gas. It would be true when b = a, and also when b ^ a if the bound systems 
are out of thermal equilibrium with free a and b, however it is not always the case 1121. We do not see the 
opposite since it is not of principle here. 

One can define dQ through the energy lost during a6-scattering. If the pair is combined (impact param- 
eter p < /9 max ), then their energy = E a + Eb = + rn, T> j s l os t completely. Otherwise (p > p ma , x ), 

their energy is lost partially due to dipole radiation, what is for given p and v = \v a — Vb\ (in large scattering 
angle limit, which is realized when p <C a y /(pv 2 )) 
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where E re \ = pv 2 /2 is the energy of relative motion. So, energy losses rate by ab -gas per unit volume can 
be given by 


e = n a n b 


rpmax f"D 'j 
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3 aa-, 66-scatterings do not lead to dipole radiation 
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where f a} b is the distribution in velocity (Maxwell). Upper limit D should be given by Debye length of 
a6-plasma, but thanks to fast convergence of AE(v,p ) with p —> oo, we put D -» oo in Eq. (18). Averaging 
over velocity distributions gives 


( ApairUrec^) — _ &0 ( m ) -^oj 


P-i _ 
2 2 o/tT 


T n 


JS-1 


(7-/3)T(2-§) 


p-i 


{EM) Ta . 

2 2 y^ \J-a J 

Combining dQ = —eV dt with other terms of the 1st law of thermodynamics gives 


3 dT a 
2 dt 


3 1 1 
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( 21 ) 


Note, that the first two terms in the right side of Eq.(21), oc — ji^pair oc (/? — 1), originate from 
the fact of disappearance (recombination) of the pair and lead to a heating of ab -gas (at (3 > 1). Term 
oc ^Eei does not dominate for a rec given by Eq. (fT|) and diminishes this effect. It is clear after accounting 
for Eqs.(fl6l).(fT9l),([20l) 
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Then Eq.(21) reduces to the second equation of the system ([6]). Hubble parameter at RD-stage is H = 

\J w ith m Pi being the Plank mass. At MD-stage it can be roughly given by H(RD) i/Trm/T , 
where Trm 1 eV is the temperature when RD—>MD transition occurs. The late A-dominated stage is not 
considered. 

The system Q can be transformed to 


dr r 2 

dT s 


31 = -ry£) 


r9 2 
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dT ’ s fps+1 ' 

Here the following notations are introduced: index s =’R’ or ’M’ means RD- or MD- stage, 

7= 5 ^ 11 , /3r = /3 - 2, /3 m = /3-|, 


Dr = Dm\/Trm = 


r (2 - f) 9s 

2 2 Y^9e 
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a 0 mpi ( pT ay ) 2 


To solve the system (24), one divides second equation by first, from where one gets independently on the 
stage 

r(9) = rp (E 1 / 7 , (25) 
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for initial conditions 0(T$ = T ay ) = 1, r(To) = r$. Being interested in r(T) on MD-stage, we will put the 
solution of system for RD-stage to be initial conditions for that for MD-stage. Substituting r(9) of Eq.(25) 
in this manner into second equation of the system (24) yields 
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9(T) = 1 + ^ R , 0 
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where 7 = 1 + 7 ^ 7 ^. Function 0(T ), with /3 = 14/5, is very slowly growing with decrease of T (/3 R = 4/5, 
/3 m = 3/10, 7 = 1/6, 7 = 23/20). At T < T RM < T 0 , 0(T) oc T" 1 / 23 . 


Solution for r(T) is given by Eqs.(25) and (26). At T <C Trm To (with Pr.m > 0) 

7-!/3m 
° 7 


(27) 


Note, that a weak sensitivity of final density to its initial value (r 0 ) is obliged to thermal effects (second 
equation of (J6]) ), and it vanishes totally if to ignore them ( 0 = 1 , 7 = 0 , 7 = 1 ). 

Appendix C 

The action of the pair of mutually attracted particles looks like 

*2 

S- f 


tl 


——b — ) dt. 
2 r 


(28) 


Here v is the current velocity rather than initial one as defined above. In the region of interest we have 
P/- ~ whence v ~ One replaces dt by dr 
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a y dr 


r v 


^Jj^y (V^2 - Vn) ■ 


(29) 


If we choose r\ = R^ and r 2 = Li OS s + Rb ~ 10 2 ^Rb then condition S> 1 gives Eq.Q. Condition Rf, S> as 
(and consequently Eq.Q) can be explicitly derived if we take n = as and 772 = Rb- On this interval the 
most of energy (~ I) is lost. 
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